GLOBAL (IN TIME) SOLUTIONS TO THE 3D-NAVIER-STOKES 

EQUATIONS ON 

T. L. GILL AND W. W. ZACHARY 

Abstract. In two recent papers ( [GZl) |GZ2) '). we provided solutions to the 
well-known unsolved problem of constructing sufficiency classes of functions 
in IHI[R'^]'^ and V[IR'^]'^, which would allow global, in time, strong solutions to 
the three-dimensional Navier-Stokes equations. These equations describe the 
time evolution of the fluid velocity and pressure of an incompressible viscous 
homogeneous Newtonian fluid in terms of a given initial velocity and given 
external body forces. In both previous papers, our solution was restricted to 
functions defined on a bounded open domain of class C"^ contained in R^. In 
this paper, we study this problem for functions defined on all of R'^. We prove 
that, under appropriate conditions, there exists a positive constant a and a 
number u+, depending only on the domain, the viscosity, the body forces and 
the eigenvalues of the "Hermite" Stokes operator (defined below) such that, 
for all functions in a dense set D contained in the closed ball B(R'^) of radius 
(l/2)u+ in IHI[R^]^, the Navier-Stokes equations have unique strong solutions 
in Ci ((0,oo),H[R3]3). 

Introduction 

Let L'^[M'']'^ be the real Hilbert space of square integrable functions on R'^ 
with values in M.^, and let IHIo[IR.'^]'^ be the completion of the set of functions in 
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{u G Cg°[M'^]'^ I V • u = 0} which vanish at infinity with respect to the inner prod- 
uct of L^[M'^]'^, and let Vo[M'^]"^ be the completion of the above functions which van- 
ish at infinity with respect to the inner product of IHIq[IR'^], the functions in EIIo[K'^]'^ 
with weak derivatives in (L^[M'^])'^. The global in time classical Navier-Stokes 
initial-value problem (on and all T > 0) is to find functions u : [0, T] x ^ 
and p : [0, T] X R3 _^ ]R guc^ that 

dtn + (u • V)u - i^Au + Vp = i{t) in (0, T) x 
V • u = in (0, T) X R3 (in the weak sense), 

(1) 

lim u{t, x) = on (0, T)xR^, 

||x|| — i-OO 

u(0,x) = Uo(x) in 

The equations describe the time evolution of the fluid velocity u(x, t) and the 
pressure p of an incompressible viscous homogeneous Newtonian fluid with constant 
viscosity coefficient v in terms of a given initial velocity Uo(x) and given external 
body forces f(x, t). (Note that our third condition, lim u(i,x) = on (0,T) x 

|x|| — >oc 

R^, is natural in this case since it is well-known that H§[R''']''' — H'^[R^]^ (see Stein 

m or m-) 

Purpose 

Let P be the (Leray) orthogonal projection of (L^[R^])^ onto Ho[R^]^ and define 
the Stokes operator by: Au =: -PAu, for u e D(A) C IHg[R3]3, the domain of A. 
Let Bu =: l/2P(-A+|xp)u for u e D{B). We caU B the Hermite-Stokes operator. 
The purpose of this paper is to prove that there exists a number u+, depending 
only on A, B, /, 1/ and R^, such that, for aU functions in ID) L'(A) nB(R^), where 
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B(M'^) is the closed ball of radius u+ in Ho(M'^)'^, the Navier-Stokes equations have 
unique strong solutions in u e Li^J[0, cx)); Vo(M3)3] n Ci[(0, oo); HolM^)^. 



Preliminaries 

In terms of notation and convention, we follow Sell and You [SY] , In or- 
der to simphfy notation, we let H denote ]HIo[M^]^ and V denote Vo[M^]^. Our 
use of the Fourier transform follows the definition of Rudin [RUj : ^{h) = 
p^ja^a /b3 ^^^^ l^{y)d'y , so that no factors of 2tt appear in the transform pairs. In 
order to simplify our proofs, we always assume that all functions u, v are in D{A) 
and, as in [GZ2j . we take c = max{ci}, where Ci is one of the nine positive constants 
that appear on pages 363-367 in 'SYJ. It will also be convenient to use the fact that 
the norms of V and V^^ are equivalent in their respective graph norms relative to 



The Stokes Operator 

It is known that A is a nonncgative linear operator which generates an analytic 
contraction semigroup. It follows that the fractional powers A^/^ and A^^/^ are 
well defined. Moreover, it is also known (cf., [SY| . |Tlj ) that the norms ||A^/^u||j^ 
and II A~^/^u||jj are equivalent to the corresponding norms induced by the Sobolev 
space (i7i[M3])3^ go that: 



(2) 



Ai/2u 



and llul 



A-i/2u 



In addition, A is an isomorphism from -D(A) °"*°> D{A^^). Furthermore, the 
embeddings V ^ H ^ V^^ are continuous, and it is easy to see that A^^ is the 
projection of an operator represented by the Riesz potential, mapping D(A^^) 
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onto D{A) (sec Stein ^). Applying the Leray projection to equation (1), with 
C(u, u) = P(u • V)u, we can recast equation (1) in the standard form: 

dtu = -i^Au - C(u, u) + Pi{t) in (0, T) x R^, 



(3) 



V • u = in (0,r) X M^, 
Hm u{t,x) = on (0,T) x 

[|x|Hoc 

u(0,x) = Uo(x) in 

where we have used the fact that the orthogonal complement of ]HI[R'^] relative to 
(L^)[R'^])'^ is {v : V = Vg, q G (iJ^ [R'^])'^} to eliminate the pressure term (see 
Galdi [GA] or [SY, Tl, T2]). Theorem 1 below will be used to get our basic estimate 
in Theorem 3. This result is a simple extension of the bounded domain case first 
proved by Constantin and Foias [CFj . 



Theorem 1. Let ai,l < i < 3, satisfy < ai < 3, < 012 < 2, < as < 3, with 
ai + a2 + > 3/2 and 

(ai, a2, as) ^ {(3/2, 0, 0), (0, 3/2, 0), (0, 0, 3/2)} . 

Then there is a positive constant c = c{ai) such that 



|(c(u,v),w)jj| 



A"i/2u 



^(l+Q2)/2^ 



We shall make use of the following interpolation inequality: (see Sell and You 
SY] , page 363) 



|A^u|U^c||A"uC||A'5u||^'-''^ 



for aU u e -D(A"), where 7 = 6la + (1 - e)l3, a, /3, 7 £ R, < 6* < 1 and /3 ^ a. 
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The Hermite-Stokes Operator 

The operator B = l/2(— A + |xp) is the three-dimensional version of the stan- 
dard harmonic oscillator operator, which generates the Hermite functions (products 
of the Hermite polynomials by e^^ as eigenfunctions for the eigenvalue problem 
on R, ( see Hermite |HR] . Appell and Kame de Feriet |AK| . and Magnus, Oberhet- 
tinger and Soni |MOS| ). It is easy to show directly, by separation of variables, that 
the solution to the 3-dimensional problem is the product of the solutions to the 
1-dimensional problem, while the eigenvalues for the 3-dimensional Hermite poly- 
nomials are the sums of those for the 1-dimensional polynomials. Furthermore, B, 
and hence B = PB, is positive with a compact inverse, while A has an unbounded 
inverse on ]HIo(R'^)'^. It turns out that B is "natural" for M.'^ in the sense that it 
is the only positive self-adjoint (sectorial) operator of lowest degree that is invari- 
ant under both rotations and Fourier transformations. (This is actually true for 
M", n > 1.) 

We will have need of the fact that every function h{t) £ H has an ex- 
pansion in terms of the eigenfunctions of B so that, for example, H h.{t) = 
J2'kLi^k^^k{t)e^{x) and, from here, it is easy to see that ||B~'^ h(i)||jj < 
A^^ II h(i)||jj, where Aj~^ is the largest eigenvalue of B^^. We also need the follow- 
ing result for our basic Theorem. 

Lemma 2. D{A) = D(B). 

Proof. If we define a norm on D{A) by ||u||^ = ||Au||jj, then {D{A) , || • ||^) is 
a Hilbert space. Now note that the Fourier transform ^{■) is an isometric iso- 
morphism on {D{A) , || • ||^) to (z?(P |x|') , || • ||^) , since ||Au||h = ||;?(Au)||jj = 
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P |x|^ u . It is now easy to see that D{A) = D{P |x|^). From this, it foUows that 

M 

D{A) = D(B). □ 



It follows from the above lemma that (AB)^* is bounded for S > 0. The 
following estimate is equation 61.24.1 on page 366 in Sell and You [SYJ. If we set 
ai = 1, Q!2 = 1/2, and 03 = in Theorem 1, along with the interpolation inequality, 
we get that 



(4) 



(C(u,v),w)nKc 



Ai/2u 



1Av||h||w||ji. 



Theorem 3. Let u, v, w e H, and let e > be arbitrary. Then, for S = 1/4 + 6/2, 
we have that: 

(5) ((AB)-(i+^)C(u,v),w) <cAr^^+'^||u||j,||v||j,||w||j,. 

Proof. Using the self-adjoint property of A, and integration by parts, we have 

(A-^'CK v), h)jj = (C(u, v), A-^h)^ = - (C(u, A-'^h), v)jj . 
It now follows from Theorem 1 that: 



|(A-'3C(u,v),h)jj| 
If we set f3 = I + S, ai — — 0, we have 
/A-(i+'5)c(u,v),h 



c 


A"i/2u 




^-/3+(l+a2)/2jj 










H 




H 





^ c||u||h||v||jj 



A("2-1-25)/2j^ 



With (5 = 1/4 + e/2, we get that, for the last term to reduce to ||h||jj, we can set 
a2 = 3/2 + e. It follows that the conditions of Theorem 1 are satisfied if 3/2 + e < 2. 
Thus, it suffices to assume that e < 1/2, which we will do in the rest of the paper 
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without comment. Our proof is completed by taking h = B~^w, and the fact that 
||B-'3w||^<A-'^||w|U. □ 

Example 4. // we use Theorem 1, with ai = 5/4, a2 = 1/4, and as = 0, along 
with the interpolation inequality, and the fact that ||A^/^u||jj ^ ||Au||jj we have 
that, for all u,v G D{A), 



||C(u,v)||h<c Ai/^u 

(6) 

s$ cllAulLllAvlL. 



^iiAuiir 



A^/^v 



riiAviir 



A better estimate is possible, but for our use, equation (6) will suffice. 

Definition 5. We say that the operator J(-,t) is (for each t) 

(1) 0-Dissipative if (J(u, t), u)jj < 0. 

(2) Dissipative if (J(u, t) — J(v, f), u — v)jj < 0. 

(3) Strongly dissipative if there exists an a > such that 



(J(u,t) - J(v,i),u- v)jj < -a||u- v||h 



(4) Uniformly dissipative if there exists a strictly monotone increasing function 
a(t) with a(0) ~ 0, Umt^oo ciit) ~ oo, and: 

(J(u,t) - J(v,i),u- v)j, < -a(||u- v||jj) ||u- v||jj. 

Note that, if J(-,i) is a hnear operator, definitions 1) and 2) coincide. Theorem 
6 below is essentially due to Browder [B], see Zeidler [Zl Corollary 32.27, page 868 
and Corollary 32.35, page 887 in. Vol. IIB], while Theorem 7 is from Miyadera [Ml 
p. 185, Theorem 6.20], and is a modification of the Crandall-Liggett Theorem |CL] 
(see the appendix to the first section of [CL] ) . 
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Theorem 6. Let B[M''] be a closed, bounded, convex subset of M[R^]. If J{-,t) : 
B[M^] IHI[IR^] is closed and strongly dissipative for each fixed t>0 then, for each 
b e M[M.^], there is a u G M[M.^] with 3{u,t) = b (e.g., the range, Ranp{-,t)] D 
M[M.^]). 

Theorem 7. Let { A{t), t £ I — [0, oo)} be a family of operators defined on M.[R^] 
with domains D{A{t)) = D, independent oft. We assume that P = D nB[]R^] is a 
closed convex set (in an appropriate topology): 

(1) The operator A{t) is the generator of a contraction semigroup for each 
tel. 

(2) The function A{t)u is continuous in both variables on / x D. 

Then, for every uq € U, the problem dtu{t,x) = A{t)u{t,x), u(0, x) = uo(x), has 
a unique solution u{t, x) e (J; D) . 

M-DissiPATivE Conditions 

Let us assume that f{t) G i°°[[0,oo);H] and is Lipschitz continuous in t, with 
\\f{t) - f(T)||H < d\t-Tf , d > 0, < e < 1. With S as in Theorem 3, we can 
rewrite equation (3) in the form: 

dtu = i/(AB)i+'^J(u,t) in (0, T) x n, 

(7) 

J(u, t) = -B-(i+*)A-*u - z/-i(AB)-(i+*)C(u, u) + i/-i(AB)-(i+*)pf (t). 

Approach 

We begin with a study of the operator J{-,t), for fixed t, and seek conditions 
depending on A, B, u, and f(t) which guarantee that 3{-,t) is m-dissipative for 
each t. Clearly J(-,i) : £)[(AB)(i+*)] — £)[(AB)(i+*)] and, since z/(AB)(i+*) 
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is a closed positive (m-accretive) operator (so that — (AB)(^+''' generates a linear 
contraction semigroup), we expect that u{A'B)^^~^^^ J{-,t) will be m-dissipative for 
each t. 



Theorem 8. For t G I = [0,oo) and, for each fixed u G H, 3{u,t) is Lipschitz 
continuous, with ||J(u, — J(u, r)||][j < d'\t — Tf, where d' = du~^a~'^^~^^\ d is 
the Lipschitz constant for the function f(t) and a~^^~^^^ = (AB)"^^'*''^^ 



Proof. For fixed u € H, 



||J(u,i) - J(u,r)||H = iz-i (AB)-(i+^)[Pf(i) -Pf(T) 

<dt^-'a-^'+'Ut-Tf =d'\t-T\r 



□ 



Main Results 

Theorem 9. Let f = supj£j^+ ||Pf(f)||jj < oo, then there exists a positive constant 
u+, depending only on f. A, B and v such that, for all u with ||u||][j < u+, J{-,t) 
is strongly dissipative. 

Proof. The proof of our first assertion has two parts. First, we require that the 
nonlinear operator 3{-,t) be 0-dissipative, which gives us an upper bound u+ in 
terms of the norm (e.g., ||u||jj ^ u+ ). We then use this part, and the fact that 
||u||jj < ||Au||jj, to show that J{-,t) is strongly dissipative on the closed ball, 
B+ = {u e H : ||Au||u < (l/2)u+}. 
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Part 1) From equation (5), we consider the expression 
(J(u, t), (AB)-^u)jj = - <B-i(AB)-^u, {AB)-'u)^ 

+ v-^ (-(AB)-(i+^)C(u, u) + (AB)-(i+^)pf (t), (AB)-^u 



B-i/2(AB)-'^u 
B-i/2(AB)-''u 
It follows that 

(J(u, i), (AB)-^u)^ < - B- V2(AB)-^u 



1 ((AB)-(i+'')C(u,u),(AB)-'*u\ ((AB)-(i+'')pf(<),(AB) 



1 <JC((AB)-(i+*)u, u), (AB)-*u^^ + v-^ {^(AB)-(i+*)Pf (i), (AB) 



C((AB)-(i+*)u, u), (AB)-^u 



+ j,-ia-(i+^)/||(AB)-^um 



In the last line, we used our estimate from Theorem 3. We now choose the first 
eigenvalue A„, n > 1, and number w such that 

(1) \-^'\-^ ||u||jj < ||B-V2(AB)-*u||jj < A^'/V* ||u||u , 

(2) \--'\-' ||u|U ^ ||B-V2(AB)-^u|L < \-,"\-' ||u|| 



and let Aq ^ = max{A„^, A^ '^}. It then follows that — Aq ||u||^ > 
- ||B-V2(AB)-*u||^. Thus, J(-,i) will be 0-dissipative if 

-Ao^a-^^ ||u||^ + ca-^(.Ar^))-i ||u||^ + (.a(^+^^))- V ||u||h < 0, 



so that 



(8) a-' ||u|U UvX^^^'^r' ||u||^ - Ao^a-^ ||u|U + {va^'+'^)-' f 



^0. 



Since ||u||u > 0, we have that 3{-,t) is 0-dissipative if 



z{yX^^+'^)-' ||u||^ - Ao^a-* ||u||h + {ya^'+'^T^ f < 0. 
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Solving, we get that 



where 7 = (4cA^/) /(z.2a(i-^)A^'+^^). Since we want real distinct solutions, we must 
require that 



^ = (4cA2/) /(i/2a(i-^)Ai^+')) < 1 ^ i.2^(i-^)Ai^+') > 4cXlf 

V > 2Aoa-(i-^)/2A-(^+')/'(c/)i/2. 



It follows that, if Pf ^0, then u_ < u+ , and our requirement that J is 0-dissipative 

implies that, since our solution factors as (||u||jj — u+)(||u||jj — u_) < 0, we must 
have that: 



|u||h-u+<0, ||u||h-u_>0. 



First observe that terms of the form (AB)~*u are dense. Then note that J{u,t) 
is closed, and the dissipative nature of an operator is determined on a dense set. 
It follows that, for u_ < ||u||jj < u+, (J(u, t),u)|j < 0. (It is clear that, when 
¥f{t) = 0,u_ = 0, and u+ = z/(cAoa*)-U^^+''\) 
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Part 2): Now, for any u, v e H with max( ||Au||jj , ||Av||jj) < (l/2)u+, we have 
that 



<J(u, t) - J(v, t), (AB)-^(u - v))jj = - B-V2(AB)-^(u - v) 
- v-^ ^(AB)-(i+'') [C(u, u - v) + C(v, u - v)], (AB)-*(u - v) ] 

< -X,'a-'' ||u - v||^ + ca-'v-^X-^'+'^ ||u - v||^ (||u|U + ||v|U) 

^ \-l -25 II ||2 I -S -l\-(l+i5)ii ||2 

< — Aq a ||u — v||u + ca V ^ ||u — v||][jU+ 

= -Ao^a-^^ ||u - v||^ + ca-'.-'X-^'+'^ ||u - v||^ [\uX^^^'\c-\-'X-,') {l + 

= _iA-ia-2* ||u-v||^{l-V^} 

= -a ||u - v||h , a = ^Aq ^a"^* |l - ^1 -7} • 

□ 

Theorem 10. The operator A{t) = i/A(^+''^ J(-, is closed, uniformly dissipative 
and jointly continuous in u and t. Furthermore, for each t G R+ and /? > 0, 
Ran[I — (iA{t)\ D B[f2], so that A{t) is m-dissipative on P. 

Proof. Since J{-,t) is strongly dissipative and closed on B, it follows from Theorem 
6 that Ran[3{-,t)] D B. 

To show that A{t) = i'(AB)(^+'^) J(-, is uniformly dissipative for u, v G B_|_, 
we have 



{A{t)u - A{t)v, (u - v))h = -I. A'/\n - v) 



- ((l/2)[C(u - V, u) + C(u - V, v)], (u - v)) 
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Now, from equation (4), 



K[C(u-v,u) + C(u-v,v)],(u-v))h| 



^ c 



aV2(u - V) ||(u - v)||ij {||Au||„ + ||Av||„} . 



We now use — Aq ||(u — ^ — ||A^/^(u — v)||jg, and the fact that the first 
eigenvalue of B is 1/2, so that A^^'' < 1, to get: 



{A{t)u-A{t)v,u-v)^^-iy Ai/2(u-v) 



+ 5C 



aV2(u - v) ||(u - v)||h {||Au|U + ||Av||h} 





Ai/2(u- 






AV2(u- 


v) 




AV2(u- 


V) 



vlljiillAullH- 



lAv 



y} 



AV2(u - V) II ||u - vm {-1 + yi^} < 0. 

1 1 IHI ^ J 



If we set a(||(u- v)||h) = -^lyX^^a-^ + ||aV2(u - v) we have 

that: 

{A{t)u - A{t)v, u - v)h < -a (||(u - v) || (u - v) ||h • 



It follows that A{t) is uniformly dissipative. Since — A(^+'*) is m-dissipative, for /? > 
0, i?an(7+/3(AB)(^+''^) = H. As J is strongly dissipative (in the ball of radius 5U+) 
and closed, with Ran[J] D B, and J{-,t) : D °"*°) D, A{t) is maximal dissipative 
(in the ball of radius ^u+), and also closed, so that Ran[I — PA{t)] D B. It follows 
that A{t) is m-dissipative on B for each t G R+ (since H is a Hilbert space). To see 
that A{t)u is continuous in both variables, let Un,u e B + , ||A(u„ — u)||jj — > , 
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with tn,t € I and tn t. Then (see equation (6)) 

||^(f„)u„ - ^(i)u||H < ||^(*n)u - ^(i)u||u + ||>l(t„)u„ - >l(t„)u||]g[ 

= ||[Pf(t„) - Pf(t)]||M + |kA(u„ - U) + [CK - U>Un) + C(U,U„ - U)]||jj 

^d\tn- tf + V ||A(u„ - U)||jj + ||C(u„ - u, u„) + C(u, U„ - u)||jj 

< |t„ - tf + V ||A(u„ - u)||jj + c ||A(u„ - u)||h {||Au„||h + ||Au||jj} 

< d |i„ - tf + V II A(u„ - u)||jj + +2c II A(u„ - u) ||jj u+. 

It foUows that A{i^)\x is continuous in both variables. □ 

Since B+ is the closure of D = -D(A) fl B equipped with the restriction of the 
graph norm of A induced on £'(A), it follows that B+ is a closed, bounded, convex 
set. We now have: 

Theorem 11. For each T e R+, t G (0, T) and uq e D C B, the global in time 
Navier-Stokes initial-value problem in : 

dtu + (u • V)u - i^Au + Vp = {{t) m (0, T) x W^, 
V-u = m fO,T) xK^, 

(9) 

lim u.{t, x) = on (0, T) x K^, 

||x||->oo 

u(0,x) = uo(x) in M^, 

has a unique strong solution u(t,x), which is in ^^^^[[0, oo); H^] and in 
L^,[[0,oo);V]nCi[(0,oo);H]. 

Proof. Theorem 7 allows us to conclude that, when uo S ID), the initial value problem 
is solved and the solution u{t, x) is in C^[(0, oo);D]. Since D c tf, it follows that 
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u{t, x) is also in V, for each t > 0. It is now clear that, for any T > 0, 

r 2 2 

/ ||u(i, x)||jj < oo, and sup ||u(i, x)||y < cx3. 

Jo 0<t<T 

This gives our conclusion. □ 

Discussion 

It is known that, if Ug G V, and f (i) is L°°[(0, oo), H] then there is a time T > 
such that a weak solution with this data is uniquely determined on any subinterval 
of [0,T) (see Sell and You, page 396, |SY| ). Thus, we also have that: 

Corollary 12. For each t e R+ and Uq G D the Navier-Stokes initial-value problem 
on M? : 

dtu + (u • V)u - i^Au + Vp = {{t) in (0, T) x M^, 
V • u = in fO,r) X E^ 

(10) 

lim u(i, x) = on (13, T) x 

||x|| — oo 

u(0,x) = Uo(x) in M^. 

/las a unique weak solution u(i,x), which is in Lf^JjO, oo);M^] and in 
L-[[0,cx));V]nCi[(0,(^);H]. 

Since we require that our initial data be in H^, the conditions for the Leray-Hopf 
weak solutions are not satisfied. However, it was an open question as to whether 
these solutions developed singularities, even if Uq £ Cg° (see Giga G and references 
therein) . The above Corollary shows that it suffices that uq (x) e H'^ to insure that 
the solutions develop no singularities. 
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